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Abstract The OmegaCAM wide-field optical imager is the sole instrument on 
the VLT Survey Telescope at ESO's Paranal Observatory. The instrument, as 
well as the telescope, have been designed for surveys with very good, natural 
seeing-limited image quality over a 1 square degree field. OmegaCAM was 
commissioned in 2011 and has been observing three ESO Public Surveys in 
parallel since October 15, 2011. We use the Astro- WISE information system 
to monitor the calibration of the observatory and to produce the Kilo Degree 
Survey (KiDS). 

Here we describe the photometric monitoring procedures in Astro- WISE 
and give a first impression of OmegaCAM 's photometric behavior as a function 
of time. The long-term monitoring of the observatory goes hand in hand with 
the KiDS survey production in Astro- WISE. KiDS is observed under partially 
non-photometric conditions. Based on the first year of OmegaCAM operations 
it is expected that a ~ 1% — 2% photometric homogeneity will be achieved for 
KiDS. 

Keywords imager, wide-field system • survey system • VLT/ VST • Astro- 
WISE • information system 

1 Public Surveys of the Southern Hemisphere 

The European Southern Observatory is in the process of delivering to its com- 
munity three Public Surveys that image the Southern Hemisphere: the Kilo 
Degree Survey [3], ATLAS and VPHAS+ pQ. The surveys are performed using 
the OmegaCAM wide-field optical imager with 1 square degree FoV. It is the 
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sole instrument on the 2.6m VLT Survey Telescope at Paranal Observatory, 
Chile. 

The three ESO Public Surveys span thousands of square degrees of sky 
obtained over years of observations. For example, the Kilo Degree Survey will 
map 1500 square degrees of extragalactic sky in the Sloan u, g, r and i pass- 
bands. This requires ~ 440 "Full Night Equivalents" spread over at least 3 
years. The scientific drivers for the survey ask for ~ 1% accuracy on the pho- 
tometric scale, in absolute sense and between passbands (see [4] in this issue). 
This means that all systematic variations along the signal path (atmosphere, 
telescope, instrument) must be characterized to a fraction of that 1%, con- 
sistently over the many years of survey operation. This includes for example 
Paranal's atmospheric variability, observed to have peak-to-peak variations of 
up to ~ 0.07mag/airmass under clear sky conditions [12] ; systematic changes 
in mirror reflectivity (dusting), with typical rates of few tenths of magnitude 
per year and subtle variations in the system's electronic gain. 

OmegaCAM and VST have been commissioned in 2011 [9]. Survey oper- 
ations started October 15, 2011 with the three surveys being performed in 
parallel. Operations are based on a strict observing protocol for both cali- 
brations and surveys. We use the Astro- WISE information system to monitor 
the calibration of the observatory and to produce and deliver the Kilo De- 
gree Survey. This paper presents the approach implemented in Astro- WISE 
for photometric monitoring OmegaCAM to maximize the quality of its survey 
imaging products. In Astro- WISE, the homogeneous stream of data flowing 
from the instrument's Calibration Plan is turned into pixeldata and metadata 
through which the full physical system of Paranal's atmosphere + VST + 
OmegaCAM is continuously monitored. The design of Astro- WISE has been 
tuned explicitly to have survey production take advantage dynamically and 
very much automatically of increases in calibration accuracy and improved 
insight of the physical system from long-term trend analyses. In other words, 
the information loop between calibrating an instrument and calibrating survey 
data is closed inside the Astro- WISE system. 

The paper is ordered as follows. We start out by describing the Omega- 
CAM instrument in its commissioned form in Section [2] We give an overview 
of OmegaCAM's dataflow (Section pjj). We then present the Photometric Cali- 
bration Plan (Section |4| and the conceptual approach in Astro- WISE for pho- 
tometric monitoring based on this plan (Section [5]). This is combined with pre- 
liminary results of OmegaCAM's photometric behavior in the first ^6 months 
of operations. We end with conclusions and outlook (Section [6]) . 



2 OmegaCAM at the VST 

The VLT Survey Telescope and OmegaCAM have both been designed specifi- 
cally to have (i) the stability to produce homogeneous surveys with long-term 
observational programs and (ii) image quality as one of the major scientific 
strengths. 
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The VLT Survey Telescope [3] is of modified Ritchey-Cretien design and 
has an alt-az mount . The VST telescope can work in two configurations. In the 
standard configuration a two-lens field corrector is used. For very large zenith 
distances the second configuration replaces this corrector with one including 
an Atmospheric Dispersion Corrector (ADC), consisting of one lens and two 
counter-rotating prism pairs. The operating wavelength ranges are 320-1014nm 
and 365-1014nm for the two-lens corrector and corrector + ADC respectively. 
The VST has active primary and secondary mirrors (Figure [T]). 

OmegaCAM is mounted at the Cassegrain focus (Figure [2]) . The optical 
parts located in the instrument are the filters, and the dewar window to the 
cryostat, which doubles as a field lens (Figure [3]). In front of the dewar win- 
dow is the mechanical part of OmegaCAM: closest to the window sits the 
filter exchange mechanism, and above that the shutter. The housing provides 
the mechanical link between the telescope and the detector/cryostat system. 
OmegaCAM is a 268 Megapixel camera with a 1 square degree FoV. The opti- 
cal design of the camera has little aberration over the full field and a constant 
0.21" plate scale. Wavefront sensing is done in OmegaCAM. The wavefront 
sensors work by registering star images that are significantly out of focus: the 
two auxiliary CCDs are mounted out of the focal plane, one 2 mm above and 
one 2 mm below. The resulting PSF is observed to equal the atmospheric see- 
ing down to 0.6" over the full FoV. Paranal offers such excellent seeing much 
of the time (a 50 % fractile of 0.66" FWHlvQ. 



2.1 Calibration unit 

The VST enclosure houses an in-dome photometric calibration unit. It con- 
sists of two sets of 4 lamps each. They are commercial 12-24V halogen lamps 
mounted on the trusses on which the secondary mounting rests (Figure [I]) . 
The lamps illuminate a dome screen, that is located on the inside of the enclo- 
sure near zenith. Each set is operated independently and is using a stabilized 
current supply, one unit per set. Lamps are switched on/off with a gradual in- 
crease/decrease of the current over a timespan of 3-5 minutes. Lamp intensity 
variations are observed to be smaller than 1% over a month (see Section pH). 
The two sets are operated at different rates, one set ~ 0.75 hour/day the 
other set at ~ 0.75 hour/week. This ensures continuous calibration coverage. 
The calibration unit serves as an internal standard candle for a range of in- 
strument monitoring and calibration activities. Dome screen observations are 
used to perform a daily overall "Quick Health Check", to measure pixel-to- 
pixel sensitivity variations, to measure CCD linearity and to check shutter 
timing accuracy. 



1 http: / / www.eso.org/ sci / facilities / paranal/site/paranal.html 
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Fig. 1 The VST from above. The Ml cover is open and the primary mirror is visible with 
its baffle. On top, the back side of the hexapod driving M2 is visible. The four sets of two 
domelamps are also visible, mounted on the white trusses forming a dodecagon. 



2.2 Shutter unit 



The exposure shutter covers an aperture size of 370x292mm. It consists of 
two carbon fiber blades. They are driven by 2 micro-stepper motors and move 
smoothly. These movements are controlled such that each individual CCD pixel 
'sees' the opening edge of the one blade and the closing edge of the other blade 
with an identical time difference, even if the blades are still accelerating. This 
provides an impact-free, high-accuracy photometric shutter. Commissioning 
tests have shown that for an exposure time of 0.2sec / lsec, deviations from a 
homogeneous exposure are below 2% / 0.2% over the whole field of view. This is 
based on observed variation in the flux level ratio between domeflat exposures 
as a function of position in the detector mosaic (Figure [4]). The domeflats are 
taken under identical conditions, but with short and long exposure times. This 
gives an upper limit to photometric effect of shutter motion inhomogeneity as 
several other sources contribute to variations in flux level ratio (e.g., calibration 
unit photometric variations, non- linear behavior of detectors) . 
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Fig. 2 The VST as seen from the observing floor. The OmegaCAM instrument and its 
control cabinets are in the lower section, interfaced to the telescope at the Cassegrain focus. 



2.3 Filter unit 



The niters are stored in two magazines which can move up and down on 
either side of the focal plane, through large shafts in the housing (Figure |3|. 
A linear stage slides filters into the beam, where they are locked into place 
by means of movable notches. The total number of motors used in the filter 
exchange and positioning mechanism is seven. Commissioning tests verify that 
the precision of filter positioning is good. Variations in photometry due to 
optical imperfections in the filters are less than 0.3% over ~ 5 months of 
operations. This has been verified for the u, g, r and i filter. The conclusion is 
based on inspection of ratios in flux level of domcflats observed over timescales 



of months. The long-term monitoring of domeflats is presented in Section 5.2 



and results are shown also in Figure 13 To minimize operational overhead time 
the filter exchange unit is built in such a way that it allows one filter to be 
pulled into the beam while the previous one is pushed out. The filters are large 
and heavy: when fully loaded with 12 filters, the instrument contains 40 kg 
of filter glass alone! During the filter exchange process about 17kg are moved 
(mass of 2 filters and the carriage). The primary filter set of OmegaCAM 
are Sloan u', g', r', i' and z' filters. In addition, there are Johnson B and V 
filters for stellar work and for cross-calibrating the photometric systems, a 



6 



Verdoes Kleijn et al. 



1.5 Meter Instrument Attachment Flange 



Filter Storage 



Magazine-A 



HI 



r — \ lks* 

Filter Optics g^gj^ 



Filter Storage 




Detector Head 



Cooling 
System 




LN 2 

Reservoir 



Fig. 3 Schematic overview of the main components of OmegaCAM. The 1.5m diameter 
housing structure contains various labeled components. From top to bottom: the shutter 
followed by the filter housing. The filters are supplied by two magazines. Each magazine can 
be filled with maximally 6 filters. The large, curved dewar window (the final optical element 
in the VST design) is followed by the detectors. The dewar is cooled using a 40-1 Nitrogen 
cryostat. 



Stromgren v filter, an Ha filter consisting of 4 segments with redshifts of up to 
lOOOOkm/sec, and a segmented ugri filter for efficient photometric monitoring 
of the sky. Filters have been manufactured by SAGEM in Paris, and by Barr 
Associates in Massachusetts, and consist of 3-layer sandwiches of coloured or 
coated glass plates. The Sloan filters are interference filters. Figure [5] shows 
transmission curves for the Sloan filters. 



2.4 CCD detectors 



The OmegaCAM optical CCD imaging camera consists of a 'science array' of 
32 2k x 4k E2V 44-82 devices for a total area of 16384 x 16384 15/x pixels, 
corresponding to 26x26cm (Figure 0. The CCDs are thinned and low-noisc 
(5e— read-out noise). They are 3-edge buttable and cover the lxl square degree 
field of view of the VST with small gaps (25"-85") and have a pixelsize of 0.21 
arcsec/pixel. Around this science array lie four 'auxiliary CCDs', of the same 
format. The mosaic has been built at ESO's headquarters in Garching. The 
whole detector system is shown in Figure [6j Two of the auxiliary CCDs, on 
opposite sides of the science array, are used for auto-guiding using a guide 
star. The field is so large that also field rotation is auto-guided. The other two 
CCDs are used for wavefront sensing using real-time image analysis. 
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Fig. 4 Ratio image of flux in a 0.1 sec domeflat exposure in Sloan z divided by the flux in 
a 3.5 sec exposure. The shutter blades move along the x-axis in the same direction for both 
exposures. For each CCD the average flux per column is computed. The results for the top 
row of 8 detectors is plotted at its right location. The subsequent rows below it are offset 
by 0.05. The two dashed blue lines for each row indicate ±2% deviation. The systematics 
in the plots, discrete jumps from one CCD to another, are dominated by the non-linearity 
of the detectors. 



3 OmegaCAM Data flow from Paranal to Astro- WISE 

The daily output of OmegaCAM is currently ~ 70Gbyte of raw imaging data 
on average, corresponding to 140 exposures. This includes day and night-time 
calibration and science exposures. This data stream of over 2 Terabyte per 
month is automatically pulled into the Astro- WISE data handling system, 
situated approximately 11250 kilometers from the imager. OmegaCAM data 
is transferred by ESO from Paranal Observatory, Chile to headquarters in 
Garching, Germany using EVALSO 6 a . From there data is pulled automat- 
ically to the Astro- WISE node in Groningen via monitoring (using wget in 
mirror mode). The nature of the data (bias, science, twilight, etc..) is deter- 
mined from header inspection of FITS format data in Groningen. This also 
intercepts data not to be ingested in Astro- WISE (e.g., focus runs). Pixel data 
is then decompressed. The new raw data triggers ingestion into Astro- WISE: 
pixel data is stored on dataservers and metadata (i.e., non-pixel data) in the 
database (see |17j for details) and an initial rough quality assessment is auto- 
matically performed. Via this retrieval approach we have succeeded in having 
OmegaCAM data available for processing in Astro- WISE as fast as 1 hour af- 
ter its observation. A week's worth of calibration data is processed and quality 
assessed by a single calibration scientist in about 1 day. For the KiDS science 
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Fig. 5 The throughput of the entire VST system for the broadband Sloan filters in Omega- 
CAM is shown (lower) with the passbands of the standard Johnson-Cousins UBVRI filters 
shown for comparison (upper). 



data, a Quick Look pipeline monitors for new raw data and automatically 
produces coadded calibrated images. Latest instead of default calibrations are 
automatically used if available. At the Astro- WISE node in Groningen we have 
achieved production of up to 70 square degrees of coadded calibrated science 
images in 24 hours. 



Swift transport and ingestion of new data in Astro- WISE is important for 
instrument monitoring and survey efficiency. Every week of delay in notic- 
ing degradation of nightly data quality and/or instrument performance corre- 
sponds to a loss of ~ 2% for a survey per year. The science yield and rate is 
also enlarged by fast production and analysis of end-products. For the Omega- 
CAM KiDS survey for example, one can think of swift follow-up observations 
(instead of 1 year waits for return of the same RA range) of highly competitive 
fields such as high redshifts QSOs and extremely cool brown dwarfs. 
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Fig. 6 The detector mosaic at the heart of OmegaCAM in its dewar ready for installation 
in the instrument. In addition to the 'science array' of 32 CCDs, the 4 'auxiliary' CCDs for 
guiding and image quality analysis are seen. 

4 OmegaCAM Photometric Calibration Plan 

The Photometric Calibration Plan of OmegaCA]V|^] is designed to provide a 
continuous photometric characterization of the survey system. A schematic 
overview of OmegaCAM's photometric Calibration Plan is shown in Fig- 
ure [8] The types of observations allow to establish the photometric scale of 
the OmegaCAM+VST system and atmosphere separately and the observa- 
tional cadences should allow continuous monitoring, i.e., sample the expected 
timescales of variations. The system is monitored using the calibration unit 
with a daily (Quick check) and weekly (Dome Flat-Field) frequency. These 
observations are the input to determine the stability and systematic photo- 
metric behavior of OmegaCAM+VST system. On-sky observations monitor 
the system plus atmosphere, using fixed fields with frequencies of three-per- 
night (Monitoring), nightly/weekly (Z eropoint, Twilight Flatfield) and monthly 
or longer (User-^-Key, Ilium. Com.). Thus the photometric behavior of the at- 



2 See http://www.eso.org/sci/facilitics/paranal/instrumcnts/omegacam/doc for Calibra- 
tion Plan and User Requirements. 
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Fig. 7 Schematic outline of the 32 CCDs in the OmegaCAM detector mosaic used for 
science observations. The chip numbering as used by Astro-WISE and ESO is indicated 
together with the on-sky orientation. Each detector is read out with a pre- and over-scan 
area in X and a overscan area in Y. 



mosphere as a function of time can be deduced from the combination of on-sky 
and in-dome observations. 



5 Monitoring OmegaCAM with Astro-WISE 



We use the Astro-WISE information system to handle the observations from 
the OmegaCAM Photometric Calibration Plan. A major design goal of Astro- 
WISE has been to have Public Survey production take advantage dynamically 
and very much automatically of improved calibration accuracy and improved 
insight of the physical system inferred from long-term monitoring using ob- 
servations from the Photometric Calibration Plan. In other words, we have 
closed the information loop between calibrating an instrument and calibrating 
a survey inside a single system of information. We refer the reader to [5] (this 
volume) and [TFj for a comprehensive description of this Astro-WISE design. 
The approach taken in Astro-WISE is to map the calibrational requirements 
(to be achieved by the Calibration Plan) to a data model. The data model 
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Fig. 8 Graphical table overview how the photometric calibration of both OmegaCAM+VST 
and the atmosphere is monitored. From left to right, the table indicates the requirement 
number and name as used in Astro-WISE an the Calibration Plan, the used filters (yi 
indicating user bands), the used fields (SA and polar fields). The stars indicate at which 
frequency the measurement is done. For further details, see text. 



is a map that describes for each calibration requirement which information it 
needs as input to be tested. This information is expressed as products from 
types of observations of the Calibration Plan. In other words, the data model 
maps which type of observations are needed to construct the desired informa- 
tion. This mapping is a cascade from final products back to the initial raw 
observations taken for the Calibration Plan (see Figure [9]). The data model 
is conceptual: it does not fix the content of a method, e.g., the algorithm, to 
make each step in this cascade. 

This conceptual data model is implemented as an object model. This means 
that each photometric requirement has a corresponding class or classes (e.g., 
Zeropoints — > PhotometricParameters, Monitoring — > AtmosphericExtinctionCoef f icients). 
The classes aggregate the information required to test the requirement. Ob- 
jects are specific instantiations of classes. For example, each set of exposures 
that represent a component of the Calibration Plan (i.e., Observing Blocks 
in ESO speak) is represented as one object. Objects can contain both meta- 
data and (processed) pixel data. Meta-data is defined as all non-pixel data: 
e.g., observational metadata, processing configurations, image statistics and 
quality assessment parameters. The objects have links to their dependencies 
recursively, where dependencies are other objects on which it depends for infor- 
mation input. Thus, objects representing final survey products can be traced 
back to their set of raw data dependencies from which they are ultimately de- 
rived. This is the concept that eventually links requirements on calibrational 
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Fig. 9 Overview of classes of data objects. Data objects not only contain the survey prod- 
ucts denoted by familiar names in wide-field imaging. They also carry the information how 
they, as requested target, can be created out of other survey objects, illustrated by the 
arrows. Underlying is an object model that captures the relationship between requested 
information and the physics of the atmosphere-to-detector observational system. 



accuracy in the Calibration Plan to attributes of objects stored in the Astro- 
WISE database. Trend plots of these attributes as function of time, any other 
attribute of the object, or any attribute of its dependencies is the basis of 
quality assessment and control. The Python Language is used to implement 
the object modef^] 

In essence, Astro- WISE is an integrated information accumulator for sur- 
vey data handling. Not a pipeline feeding a separate archive feeding a separate 
set of data analysis tools. It keeps the cycle processings archiving — > in-depth 
monitoring analysis around the pool of data (observational pixels and meta- 
data) in a single system. Analysis by survey team members and calibration 
scientists on the data results in an improved understanding of the instrument 
and how it should be calibrated. Changing methods is one way in which this 
improved knowledge can be captured by Astro- WISE. This does not alter the 
object model. By applying the improved knowledge (i.e., reprocessing) the 
next cycle starts to increase the level of detail to which the instrument plus 
the processing of its information is understood. 

The next sections describe in more detail the observations for each com- 
ponent in the Calibration Plan and how they are handled in Astro- WISE. 



3 The Astro- WISE Python documentation server gives at http://doc. astro- 
wise. org/astro. main. html an IT-technical description of the Astro-WISE object model. 
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Preliminary monitoring results are presented. They are based on the first ~ 6 
months of operations. 



5.1 Quick check 



40000 

5 
a 

r 35000 

<u 
> 

V 



30000 



25000 



20000 u 



t r 

#81 



n 1 r 



«. •. 



_l I I L 



40000 



35000 



30000 



25000 



20000 



- — i — I — ! — 1 — 1 — : — ' — I — ~ 

#81 



_l I l J J l i - i 



Fig. 10 Example of trend analysis based on the Quick check observations for CCDs: #81. 
The Quick check procedure observes the dome screen through the Sloan r filter in a closed 
dome during daytime. Shown are the raw detector values (dots, in ADU) with their 3a 
variation (dotted lines) around the clipped mean (dashed line). Both the lamp set in primary 
use (blue) and secondary use (red) are shown. Left: zoom in on April 2012. Right: long-term 
trend for over 10 months of operations. 



Procedure: The Quick check provides a daily check on the overall system 
health in terms of responsivity. It observes the dome screen through the 
Sloan r filter in a closed dome during daytime for 0.5 seconds. For the op- 
erations at Paranal Observatory this measurement leads to the go-ahead/non- 
conformance verdict in the day report. The Quick check procedure is a comple- 
ment to the general domeflats, which are observed weekly for all filters. The 
QuickCheckFram^] is the object model implementation of the Quick check 
component in the Calibration Plan. It stores data and metadata used in Fig- 



ure 10 and provides links to (meta-)data of its dependencies. 

Usage / preliminary results: The Quick check results provide a constraint 

on the 'in-dome' absolute photometric scale of the system: OmegaCAM, VST 



see the Python Documentation server at http: //doc. astro- 
wise. org/astro. main. QuickCheckFrame.html for the description of the QuickCheckFrame 
class. 
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Fig. 11 Scatter in absolute mean exposure level of domefiats over 3 months for each CCD 
separately. The large symbols at the end gives the observed scatter for the MOSAIC as 
a whole. These measurements imply that the in-dome photometric scale is constant to a 
level of ~l-3mmag level for most CCDs for this period. For the full ~10 months of science 
operations values have been measured in the range of ~l-6mmag for timescales of a single 
month. 
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Fig. 12 Example of a gain monitoring plot for human inspection zooming in on three 
CCDs: #81, #82 and #83. The measured gain (e~ /ADU, dots) is plotted as a function of 
8 months of time. Blue/red indicates measurements deemed to be valid/invalid. The red 
dashed lines indicate ±0.5% of the median value. The stability observed for #81 and #83 is 
typical for the whole mosaic. The variation in gain measurements for #82 has disappeared 
since an intervention on June 2, 2012. The gain procedure uses observations of the dome 
screen through the Sloan r filter in a closed dome during daytime. They are done at a range 
of exposure times spanning the full range of electron well depth. The gain is derived from a 
fit to the Photon Transfer Curve built from those observations. 



and calibration unit. Figure [10] shows the trend of raw detector levels over 
one and ten months for one CCD. It implies that the photometric scale of 
the calibration unit varies almost linear on timescale of a month with a more 
complex behavior on timescale of several months. The repeatability of the 
measurement on a monthly scale corresponds to ~ 2-5 mmag per CCD on the 
system photometric scale for most CCDs. Similarity across chip and filters of 
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residuals to the mean exposure level suggests variability of the lamp is the 
dominatant source of scatter for most chips. Two CCDs show a larger scatter 
(see Figure 11 ). These two deviations must originate within electronics specific 
to these detectors. Previously, CCD#82 showed a larger variation. This was 



probably due to gain variations (see Figure 12). This has been solved with a 
replacement of a videoboard on June 2, 2012. The origin of the initial decline 
and subsequent flattening of the responsivity over ~ 10 months is not yet fully 
understood. Similar behavior has been observed for VLT UT1 and FORS2. A 
possible reason for this is that the Ml coatings degrade and the Ml becomes 
dusty (the cause of the primary decline), but after a few months this reaches 
a kind of saturation level by which the Ml cannot become much dirtier (ESO 
private comm.). 



5.2 Dome flat-fielding 
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Fig. 13 Left: mean master dome Hatfield image through Sloan i from the weekly observa- 
tions over the period 1 August 2011 - 1 February 2012. Flatficlds are per chip normalized 
to their median. The greyscaling goes from 0.98 (black) to 1.02 (white). Right: standard 
deviation image from the same set of domeflats used left. The grey-scaling goes from 1.0e-3 
(black) to 1.6e-3 (white). The typical standard deviation corresponds to < 2 mmag built-up 
out of systematic and variable changes in the system excluding the lamp. Two chips (#82 
and #92) show enhanced variation. Furthermore variations due to changes in bad columns 
and moving dust can be observed. 




Procedure: Dome flatficlds are being measured for the 4 keybands and the 
key-composite filter at least once per week. Thus at least within 3 days of the 
taking of science data a dome flatfield in the key passband is available. The 
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DomeFlatFram^jis the object model implementation of the Dome flat-fielding 
component in the Calibration Plan. 

Usage / Preliminary results: The sequence of master dome flatfields in 
the keybands, acquired over the months and eventually years is being used 
to monitor the long term stability of the instrument and calibration unit (see 



Figure 13). From the first 6 months of operations the standard variation in a 
pixel value of the master domeflat is typically 1.0e-3 - 1.6e-3 in relative units. 
Changes in bad columns and moving dust are observed to cause standard de- 
viations up to ~3e-3. In conclusion, the variation in master domeflats over 
6 months of operations suggest pixel-to-pixel variations in response to corre- 
spond to l-1.5mmag typically on the photometric scale. The variations also 
put a constraint on the per pixel photometric scale variability of the system 
without the calibration unit, as the latter one is divided out by normalization. 



5.3 Twilight flat-fielding 

Procedure: Like domeflats, twilight flatfields are being measured for the 4 
keybands and the key-composite filter at least once per week. Exposure times 
vary to yield roughly constant exposure levels, i.e., compensating twilight 
brightness variation. The TwilightFlatFram^Jis the object model implemen- 
tation of the Twilight flat-fielding component in the Calibration Plan. 
Usage / preliminary results: The master twilights acquired over the first 
~ 6 months of operations show variations higher by a factor 2 to 3 compared 



to the domeflats (see Figure 14 1. Twilight flats are observed currently at fixed 
positional angle on the sky and thus at different angles of rotation between 
instrument and telescope. A significant contribution to the variation observed 
in flatfields as a function of time is made by systematic changes in the flatfield 



as a function of rotator angle (see Figure 15 ). A possible cause is a variation in 
gradients in the stray-light distribution as a function of angle between Omega- 
CAM and VSlQ This forms the major contribution to the observed standard 



variations shown in Figure 14 . 



5.4 Zeropoints 

Procedure: Zeropoints of the overall detector chain, separately for each CCD 
chip, in at least all four keybands are obtained nightly irrespective of the 

5 see the Python Documentation server at http: //doc. astro- 
wise. org/astro. main. DomcFlatFrame.html for the description of the DomeFlatFrame 
class. 

6 see the Python Documentation server at http ://doc. astro- 
wise, org/astro. main. TwilightFlatFrame.html for the description of the TwilightFlatFrame 
class. 

7 See movie of residuals between twilight flatfield at a given rota- 
tor angle and the average of an ensemble average at http://wiki. astro- 
wise. org/projects:omegacam:portal:illuminationcorrection . 
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Fig. 14 Left: standard deviation image for the ensemble of master domeflats in Sloan g 
observed over the period 1 August 2011 - 1 February 2012. The grey-scaling goes from 1.0e-3 
(black) to 1.6e-3 (white). The typical standard deviation corresponds to 0.0013 ± 0.0003 for 
all CCDs except #82 and #92 which have a mean of around 0.0015. Furthermore variations 
due to changes in bad columns and moving dust are detectable. Right: standard deviation 
image from the emsemblc of twilightflats observed over the same period as the left image. 
The grey-scaling goes from 2.0e-3 (black) to 6.0e-3 (white). The typical standard deviation 
is two to three times higher than for domeflats. It is 0.003 ± 0.001 for the CCDs on the 
mosaic border and goes up to 0.004 ± 0.002 for the CCDs in the middle. Variations due to 
changes in bad columns and moving dust can be observed barely. An additional difference 
with the domeflats behavior are the higher standard deviations observed near the top and 
bottom gaps of the three horizontal gaps between the four rows of CCDs. This might be 
due to scattering and/or vignetting by imperfectly centered masks over the CCD connectors 
located in these two out of the three gaps (see Figure ml . 



science programme. The pointing is selected from a set of 8 equatorial SA 
fields that cover the full RA cycle (see Table [T]). Currently we use as standard 
catalogs the SDSS DR 8 data in these SA fields and the catalogs from our 
preparatory programme of Secondary Standards |16j . They will be replaced 
in the future by the OmegaCAM Secondary Standards Catalog on these 8 
SA fields. For this, a dedicated OmegaCAM observational programme is in 
progress. The PhotometricParameter^] is the object model implementation 
of the Zeropoints component in the Calibration Plan. The PhotRef Catalog 
is the object model implementation of the photometric reference catalog. 



8 see the Python Documentation server 
wise.org/astro.main.PhotometricParameters.html for the 
PhotometricParameters class. 

9 see the Python Documentation server 
wise.org/astro.main.PhotRefCatalog.html for the description 
class. 



at http : //doc . astro- 
description of the 

at http ://doc . astro- 
of the PhotRef Catalog 
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Fig. 15 Difference image between a single master twilight flat and the average of an en- 
semble of twilight flats observed at various rotator angles. The grcyscaling goes from -0.002 
(black) to +0.002 (white). Sequences ordered by rotator angle show systematic behavior 
of the differences with this angle (see text for reference to movie). This forms the major 
contribution to the observed standard variations shown in Figure [14] 

Usage / preliminary results: We model the photometric behavior of Omega- 
CAM with the common photometric equation for astronomical imagers: 



m inst = -2.5 \og{countrate) + ZPT - k x AM (1) 

where mj nst is the magnitude of the object in the instrumental photometric 
system, the countrate is in ADU/^J , k is the atmospheric extinction coeffi- 
cient, AM is the airmass and ZPT the zeropoint. We neglect color terms at 
the moment. Table [2] lists the observed scatter in on-sky measurements of the 
zeropoint, ZPT on timescales of an hour. This is inferred from series of back- 
to-back SA field observations with chip-sized offsets. This gives an upper-limit 
on how well we can measure the photometric scale. If the response of VST + 
OmegaCAM + atmosphere did not vary over the observations, this equals the 
standard deviation in the photometric scale. Thus it also provides a limit on 
the contribution of systematic error on the photometric scale. For example, we 



To convert to e /s: the typical gain for OmegaCAM;s detectors is ~ 2.5-2.6e /ADU. 
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Table 1 Standard star fields for zeropoint determination. The number of Landolt [10] and 
Stetson 14 standard stars within the OmegaCAM FoV are listed. In addition the numbers 
are listed for available SDSS (DR7) stars and stars from our Preliminary Program (PP) for 
secondary standards. 



Field 


RA [dcg] 


Dec [dcg] 


#Landolt 


#Stetson 


#SDSS 


#PP 


SA92 


13.946 


+0.949 


41 


213 


1094 


6475 


SA95 


58.500 


+0.000 


45 


426 


1093 





SA98 


103.021 


-0.328 


46 


1116 





23840 


SA101 


149.112 


-0.386 


35 


117 


1776 


5591 


SA104 


190.488 


-0.529 


34 


76 


1576 


5701 


SA107 


234.825 


-0.263 


28 


728 


3889 


12006 


SA110 


280.600 


+0.346 


39 


589 





38562 


SA113 


325.375 


+0.499 


42 


483 


4046 


13947 



assume default atmospheric extinction coefficients in this particular case. Fur- 
ther contributions can come from systematics in the reference catalog (SDSS 
was used here) and the applied illumination correction. The hourly scatter in 
the on-sky photometric scale per CCD is ~2 -7 times larger than the monthly 
scatter observed for the in-dome photometric scale. A lower-limit to how well 
we can measure the photometric scale is given by the zeropoint measurement 
precision. Table [2] shows that the current precision of the on-sky measurement 
is 1-3 times larger than the monthly in-dome scatter. The on-sky measurement 
can be improved by improving the standard catalog (SDSS DR8 catalogs in 
the SA fields at the moment) and improving the source extraction (aperture 
magnitudes with fixed 6.4 arcsec diameter, regardless of seeing). In conclusion, 
the best ZPT per-chip at any time can be obtained using in-dome observations 
with the calibration unit. However, this scale is not tied to any photometric 
system. This might be achieved by scaling the in-dome zeropoint ZPT-D us- 
ing a set of in-dome and on-sky measurements obtained over many months as 
constraint. In figures 16 and 20 an illustrative scaling (not a fit!) is overplot- 
ted on the on-sky ZPTs. It gives hope that contamination from variability in 
atmosphere and lamp can be overcome in tying the in-dome photometric scale 
to a photometric system in the near future. 



5.5 Atmospheric monitoring 

Procedure: the OmegaCAM Calibration Plan includes Atmospheric Monitoring 
of a field near the Equatorial South Pole. Every night, OmegaCAM observes 
three times the same field near the equatorial South pole: at the beginning, 
middle and end of a night. This is done through a composite filter containing a 
Sloan u, g, r and i quadrant. Table [3] lists the observational specs for the polar 
field. Additionally a polar standards photometric catalog is made from dedi- 
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Fig. 16 Photometric scale in Sloan r on sky and in dome as a function of over 10 months 
of survey operations for CCD #81. The green filled circles denote zeropoints obtained from 
SA field observations. The blue and red filled circles are dome-screen observations from 
Quick check from primary and secondary set respectively. The latter have their average 
value converted to an arbitrary zeropoint (not a fit!). The similar decrease in zeropoint 
in time indicated by all three measurements suggests a common cause. See text for more 
information. The corresponding plot for all 32 CCDs is given AS appendix (Figure [20| 



cated observations. The PhotometricExtinctionReportp^jand PhotSkyBrightnes 

are the object model implementations of the Monitoring component in the 
Calibration Plan. 

Usage: These observations serve two main goals. First, they are used to mea- 
sure the stability of the nightly atmospheric extinction. This is done by combin- 
ing information from the Zeropoints observations on SA fields and assuming 
standard atmospheric extinction coefficients. Thus it is possible to determine 
3/night the actual atmospheric extinction compared to photometric condi- 



tions (Figure 17). Second, the observations are used to measure sky brightness. 
Again these can be compared to expected values under photometric conditions 
for given lunar phase/distance. Long-term monitoring of both the atmospheric 
extinction results and sky brightness results can yield deeper understanding 



11 see the Python Documentation server at http: //doc. astro- 
wise, org/astro. main. PhotometricExtinctionReport.html for the description of the 
PhotometricExtinctionReport class. 

12 see the Python Documentation server at http: //doc. astro- 
wise, org/astro. main. PhotSkyBrightness.html for the description of the PhotSkyBrightness 
class. 
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Table 2 On-sky and in-dome accuracies and precisions in measuring the zeropoint. ZPT 
is defined as in Equation^ ZPT-D denotes the in-dome measurement of a zeropoint taking 
the calibration lamp as a 'standard candle' reference. Stdev(ZPT)MOSAlC is the standard 
deviation in zcropoints for the 32 CCDs together. Stdev(ZPT)ccD i s the same measurement, 
but now per CCD. Err(ZPT)cco is the zeropoint measurement precision derived from the 
observed scatter of calibrators around the zeropoint. (The Err(ZPT)MOSAlC is v32 ~5.7 
smaller) . 



system 


Quantity 


timescale 


u 


g 


r 


i 








(mmag) 


(mmag) 


(mmag) 


(mmag) 


on-sky 


stdev(ZPT) MO SAlC 


hour 


11 


10 


12 


11 


in-dome 


stdev(ZPT-D) MO SAlC 


month 


4 


4 


3 


2 


on-sky 


stdev(ZPT) CCD 


hour 


37 


10 


16 


~15 


in-dome 


stdev(ZPT-D) CC D 


month 


6 


4 


3 


2 


on-sky 


Err(ZPT)ccD 




15 


5 


5 


7 


in-dome 


Err(ZPT-D)- CCD 




« 1 


« 1 


<< 1 


« 1 
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Fig. 17 Inspection plot for derived excess extinction (blue points) which is the excess above 
the standard atmospheric extinction (dashed line). 



in the physics setting the atmospheric extinction (dust, aerosols, water vapor 
etcetera) . Eventually, this should lead to survey calibrations without the need 
to assume default atmospheric extinction coefficients. 

As the maximum difference in airmass is 8% between any two polar field 
observations it directly measures changes in the atmospheric extinction coeffi- 
cient at fixed high airmass to this accuracy. Similarly, comparison of overlap- 
ping sources in dithered science observations measure changes in atmospheric 
extinction coefficients at fixed lower airmasses. Results from Sections 15.11 and 
|5.4| show that the zeropoint ZPT could be constrained independently by com- 
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Table 3 Polar field monitoring specifications 



Topic 


Description 


Coordinates 


a =51.429167, <5=-89.0425 (J2000) 


Filter 


4 quadrant glass filter: u,g,r, and i 


Exposure time 


100s 


Airmass 


2.3-2.5 



billing information from Quick health and Zeropoints observations. The 
combination of information from SA fields (Zeropoints) at typical airmass~ 
1.2 and polar field (Atmospheric Monitoring) at airmass~ 2.4 can then con- 
strain the behavior of the atmospheric extinction coefficient k as a function of 
time from a nice baseline in airmass. 



5.6 Illumination characterization 

Commissioning analysis has shown that twilight and dome flat exposures with 
OmegaCAM are both subjected to non-uniform illumination. Thus, the deriva- 
tion of a 'true' fiatfield (i.e., only the variation in pixel sensitivity) requires 
a correction of these illumination variations. This is termed illumination cor- 
rectiorp^] Here we give a summary. If no correction for these illumination 
variations is applied the photometric scale residual can vary by up to ~ 0.1 



mag, bottom-to-peak within a single chip (Figure 18 1. The large-scale residual 
pattern appears roughly point symmetric around the center of the mosaic. A 
variation in the stray light distribution which is point symmetric around the 
optical axis to first order can cause this. Other instruments show similar illumi- 
nation variations patterns, i.e., having a central- axisymmetric 2D-polynomial 
shape, for example, WFI QI] and [8] and MEGACAM [B]. In addition there 
are localized components, especially near the covers over the detector heads. 
Finally, commissioning analysis has shown a dependency on the rotator angle 
between the camera and the telescope. From trend analysis on ^5 months of 
weekly dome flats and twilight flats the conclusion is that the intrinsic flat 
field (i.e., pixel sensitivities) varies with a a ~ 0.2% or less. Therefore it was 
decided to build a single masterflat from a combination of a single twilight 
flat and single dome flat for the Sloan g, r, i and z filters and apply this to 
all OmegaCAM observations used to determine the illumination correction^] 
The usage of a single flat (i.e., at single rotator angle) simplifies the procedure: 
a single illumination correction needs to be derived. 



13 A detailed description of the illumination characterization and correction is given 
in commissioning document VST-TRE-OCM-23100-3608 (available at http://wiki. astro- 
wise. org/projects:omegacam:portal:illuminationcorrection) 

14 For u only the twilight flat is used see commissioning doc- 
ument VST-TRE-OCM-23100-3608 (available at http://wiki.astro- 
wise.org/projects:omegacam:portal:illuminationcorrection) for an explanation 



OmcgaCAM and Astro- WISE 



23 



To characterize the illumination variations SA fields were observed at ~ 33 
dither positions. These dither positions are mostly single-chip offsets such that 
a group of stars is observed on all 32 CCDs. Zeropoints per chip are computed 
and the magnitude residuals, i.e., the difference between the zeropoint and 
reference magnitude, sample the illumination variations. The combined SA 
fields, yield roughly 1000 useful magnitude residuals per chip for g, r, i ,z 
and about 500 for the u-band. The residuals are modelled as a continuous 
locally linear model using adaptive binning and interpolation. The results are 
summarized in Tableland shown in Figure [19| 

Table 4 Illumination Correction Results. Standard deviations of residual magnitudes 
in the five SLOAN filters left after applying the localized illumination corrections. The inter- 
nal and external standard deviation (sd) of the magnitude residuals are defined as follows: 
external is with respect to the stellar reference magnitude in the standard star catalog 
(SDSS DR7 stellar photometry). Internal is w.r.t. the mean of its observed magnitudes. 
2nd and 3rd column: sd of internal and external errors if only zero points are applied; 4th 
and 5th columns: sd of internal and external errors when illumination correction applied on 
OmegaCAM KiDS survey data (SA fields for z band). 





zp only 


zp only 


ilium. corr. 


ilium. corr. 


(1) 


(2) 


(3) 


(4) 


(5) 


band 


sd_int 


sd_ext 


sd_int 


sd_cxt 




(mmag) 


(mmag) 


(mmag) 


(mmag) 


u 


29 


36 


13 


25 


g 


20 


30 


9 


24 


r 


21 


29 


12 


22 


i 


22 


30 


10 


24 


z 


27 


33 


15 


27 



The illumination correction is characterized to better than 1% for the am- 
plitude over a single CCD. The correction method overcomes the rotator angle 
dependence by virtue of the stability of the intrinsic ("true") flatfield. This illu- 
mination correction has remained accurate on timescales of at least 7 months 
of OmegaCAM survey operations (including a VST technical intervention). 
The illumination correction approach here can account for the sub-percent 
changes in pixel-to-pixel sensitivity over time. This can be done by updating 
the domeflat because it sets the small-scale structure in the masterflat. 

The observed flatfield and the illumination correction derived for it form a 
multiplicative pair. Together they describe the true pixel sensitivity over the 
field. The multiplication procedure does not remove the additive component 
of straylight itself. In other words, the correction should be applied to sources 
that are small compared to the scale over which variation in straylight strength 
is significant. Sources with more extended surface brightness distributions on 
scales of a CCD or larger in science exposures are still contaminated by stray- 
light. In this case a combination of a global illumination correction plus local 
background fitting appears able to remove the straylight contamination [7j. If 
no such dedicated correction is applied the extended surface brightness distri- 
butions are affected systematically up to 0.04mag for typical cases. 
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Similar to flatfielding and zeropoints the illumination correction is repre- 
sented in the Astro- WISE data model. The IlluminationCorrectiorF^lis the 
object model implementation of the IlluminationCorrection component in the 
Calibration Plan. 



6 Conclusions and Outlook 

OmegaCAM at the VST has been commissioned successfully, meeting specifi- 
cations and started science operations October 15, 2011. The continuous flow 
of observations from OmegaCAM's Photometric Calibration Plan are being 
processed in Astro- WISE. Astro- WISE is also used for photometric monitor- 
ing and analyses. The first half a year of photometric calibrations yield the 
following preliminary results: 

— The on-sky measurements with OmegaCAM achieve lOmmag precision in 
sampling the photometric scale of OmegaCAM+VST+atmosphere. 

— The in-dome photometric scale of detector + calibration unit is stable over 
a month to within ~5mmag for ~90% of the detectors. 

— Master domeflat variations over 6 months of operations constrain the rel- 
ative pixel-to-pixel responses to be constant to a level of l-1.5mmag. 

— Measuring the VST+OmegaCAM response with in-dome observations (dome- 
flats) can be done with ^2-5 times higher precision than with on-sky ob- 
servations (standard fields). First results give hope that by combining the 
results from in-dome and on-sky observations we can model variations in 
the atmospheric extinctions coefficients very accurately. 

The Astro- WISE system is also being used also to handle the Kilo Degree 
Survey. Observing conditions are partially non-photometric. The goal of our 
OmegaCAM photometric calibration is to reach photometric homogeneity to 
1% accuracy per filter and between filters. KiDS observes 1500 square degrees 
of sky in two patches. The area is tiled with ~1 square degree observations that 
have some overlap. The tiling is identical for each filter. These tiles are con- 
structed per filter from an ~ 1 hour observing block of dithered observations, 
4 in u and 5 in g, r and i. This survey design yields constraints on variations 
in photometric scale in several ways. Variations in the photometric scale on 
the hour scale are sampled by dither-to-dither comparisons within an observ- 
ing block. The overlaps between KiDS tiles provide intra-night and inter-night 
constraints. The ATLAS Public Survey running in parallel on OmegaCAM at 
the VST covers KiDS survey area in the same filters to a shallower depth, using 
a different tiling. The overlaps within KiDS and between KiDS and ATLAS 
provide therefore additional mostly inter-night constraints. Exploiting these 
constraints to improve survey photometric homogenization is an on-going ef- 
fort. 



15 see the Python Documentation server at http: //doc. astro- 
wise. org/astro. main. IlluminationCorrection. html for the description of the 
IlluminationCorrection class. 
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It is the combination of the total volume of photometric information over 
years of calibrations of OmegaCAM and survey observations of KiDS and 
ATLAS that can lead to an internally consistent continuous description of 
the observatory's photometric scale to the best accuracy. The Astro-WISE 
system design is aimed specifically at combining such large volumes of complex 
information dynamically, improving the solution as science and calibration 
observations progress and new information accumulates in Astro-WISE. 
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Fig. 18 Photometric residual of stars in SA113 after subtraction of an individual zeropoint 
per chip. Top-left to bottom-right: Sloan u,g,r,i and z. Residuals are computed using mag- 
nitudes for stars listed in the SDSS DR7 catalog. Residuals are in AB magnitudes in the 
OmegaCAM instrumental system. No color correction has been applied yet to convert the 
SDSS DR7 catalog to this system. 
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Fig. 19 Illumination correction results for Sloan r projected on the focal plane. The data- 
points in the top four plots represent all DR7 reference stars of SA92, SA95 and SA113 after 
dithering 33 times and filtering (see text). Top left: raw zeropoints for all SDSS DR7 ref- 
erence stars. Top right: model = fitted zeropoints per chip + illumination variation model. 
Middle left: residual magnitudes of reference stars after applying illumination correction. 
Middl e rig ht: illumination variation model (i.e. with ZPT per chip subtracted, similar to 
Figure [l8] l . Bottom: the resulting illumination correction pixelframe as applied in the pho- 
tometric data processing. 
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Zeropoint [mag] 




Fig. 20 Photometric scale in Sloan r on sky and in dome as a function of over 10 months of 
survey operations for all 32 science CCDs. The green filled circles denote zeropoints obtained 
from SA field observations. The blue and red filled circles are dome-screen observations from 
Quick check from primary and secondary set respectively.The latter have their average value 
converted to an arbitrary zeropoint (not a fit!). See Section [5^4] for a discussion. 



